Band dispersion and electronic lifetimes in crystalline organic semiconductors 
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The consequences of several microscopic interactions on the photoemission spectra of crystalline 
organic semiconductors (OSC) are studied theoretically. It is argued that their relative roles can be 
disentangled by analyzing both their temperature and their momentum/energy dependence. Our 
analysis shows that the polaronic thermal band narrowing, that is the foundation of most theories 
of electrical transport in OSC, is inconsistent in the range of microscopic parameters appropriate 
for these materials. An alternative scenario is proposed to explain the experimental trends. 
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Introduction. It is now becoming possible to measure 
Angle-Resolved photoemission (ARPES) spectra of or- 
ganic semiconductors (OSC) with increasingly high ac- 
curacy, revealing intrinsic properties of conduction elec- 
trons in these materials. Clear energy-momentum dis- 
persions, indicative of the formation of well-defined elec- 
tronic bands, have been reported in pentacene monolay- 
ers crystalline films 0, Q as well as rubrene single 
crystals |4|. Since photoemission probes the very nature 
of electronic excitations, the measured spectra (in partic- 
ular, their deviations from non-interacting Bloch bands) 
provide direct information on the interactions affecting 
the electron motion. Ultimately, such information could 
help establishing a proper microscopic model for the elec- 
tron dynamics in these materials. For this procedure to 
be viable, however, one must be able to provide reliable 
predictions from the theoretical side, identifying charac- 
teristic signatures of the different microscopic phenomena 
at work. Only then the comparison with the experimen- 
tal ARPES spectra can be effectively used to sort out the 
relevant microscopic interactions in crystalline OSC. 

In this work, we consider the effect of several micro- 
scopic interactions that have been considered to play a 
role in crystalline OSC: (i) the coupling of electrons to 
intra-molecular vibrations; (ii) dynamical disorder origi- 
nating from thermal inter-molecular motions; (iii) static 
disorder, of both chemical and structural origin, giving 
rise to spatial variations of the molecular site energies 
and (iv) details of the electronic structure, most impor- 
tantly the features arising from the non-equivalence of 
the molecules in the unit cell. On general grounds, these 
interactions can affect either the band dispersion, the 
electronic lifetimes (i.e. respectively the position and 
width of the ARPES peaks), or both. In order to dis- 
entangle all these effects experimentally, it is important 
to characterize their temperature dependences: for exam- 
ple, (i, ii) can in principle depend on temperature, via the 
thermal changes in the phonon population while (iii,iv) 
should be essentially temperature independent, being de- 



termined by the structural characteristics of the sample. 
A similar assessment can be made concerning their mo- 
mentum/energy dependence, i.e. how different interac- 
tions affect distinct regions of the electronic spectrum. 
Our aim here is to clarify and quantify these aspects by 
performing a reliable calculation of the electronic spectral 
function in a controlled microscopic model. 

Model. We consider the following Hamiltonian, 
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describes the tight binding motion of holes in 
a HOMO band that we take to be one dimensional 
for simplicity, reflecting the marked anisotropy that is 
commonly observed in crystalline OSC. H^> and H' n > 
are respectively the "Holstein" interaction with intra- 
molecular deformations Xj of frequency fig, with cou- 
pling strength g, and the "Pcierls" or "Su-Schricffcr- 
Heeger" (SSH) interaction with molecular displacements 
Uj of frequency ujq. These two interactions modulate re- 
spectively the molecular energies and the inter-molecular 
transfer integrals. For the latter we assume a linear de- 
pendence on the inter-molecular distance, f(uj — Uj+i) = 
Jassii(uj — u j+i)i being assH the coupling strength. 
jj(m) describes disordered site energies ej obeying a 
gaussian distribution of fixed width A, and an 
alternating potential ±5 reflecting a non-equivalence of 
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the sites in the unit cell, as is commonly observed in 
OSC. H^ mb ' is the Hamiltonian for free harmonic vi- 
brations. In the following we set h = kg = 1 and ex- 
press the strength of the electron-vibration interactions 
by introducing the dimensionless parameters an = <?/^o 
and Xssh = a |sff^/(2-^ a; o)- Typical bandwidts in or- 
ganic semiconductors are W ~ 0.5eV, corresponding to 
J ~ 0.125eV in our one-dimensional model. The micro- 
scopic parameters that can be found in the literature are 
w ^ (0.05 - 0.1)J and O ~ (1 - 2)J with X SS h ^ 
0.15 - 0.3 and a%_ ~ 0.2 - 0.5, and S = (0.1 - 0.3) J, 

a = (o.i - o.3)j ul-H, US-Eli ■ 

Method. Owing to the low frequency of the inter- 
molecular vibrations, that stems from the large molec- 
ular mass, we have both ujq <C J and ujq < T m the 
relevant temperature range. In this regime the displace- 
ments Uj in can be treated as a random field obey- 
ing a gaussian distribution, causing a statistical disorder 
in the inter-molecular transfer integrals [f| . The solution 
of all the terms in Eq. ([T]) but H^> therefore reduces 
to a problem of non-interacting electrons in the presence 
of both "site-diagonal" and "off-diagonal" disorder, that 
can be efficiently treated with the method described in 
Ref. Instead, the full quantum nature of the intra- 
molecular vibrations must be retained in i/W i because 
T « Slj. To this aim we devise a method of solution 
that provides reliable results in the case fio ^ J> which 
includes the moderate vibrational frequency regime of 
interest here, fio — (1 — 2) J. As we show below, this 
treatment of the intra-molecular electron- vibration inter- 
actions can be straightforwardly implemented in the dis- 
ordered environment provided by the remaining terms in 
the Hamiltonian. 

The Green's function of the problem can be defined in 
the site representation as 
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where P(u) oc Hi cxp(— u 2 /2<r 2 ) is the statisti- 
cal distribution of local displacements u = {w/}, 
with a 2 = [2A/a;otanh(a;o/2T)]~ 1 and P<u s (e) oc 
H m exp(— e 2 n /2A 2 ) is the distribution of local energies 
e = {e m }- For a given set of inter- molecular deforma- 
tions and disorder variables the Green's function of the 
electronic problem H e i = H — H^> in the absence of the 
Holstein term can be obtained as an inversion of a tri- 
diagonal matrix in the site representation: (u, e) = 
[(w - H el )-\j with @ 
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The interaction term ffw is subsequently included at a 
local level. In practice, we consider a thermalizcd Hol- 
stein impurity at each site i, whose Green's function is 
defined as 
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with Zh = Y^=o e nn °' T [we drop the indices (u, e) in 
the following]. Here Qi^ is the local propagator obtained 
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from the inversion of the matrix Eq. ^ emu ^ H i 

Yf^ 1 ^ are the local emission/absorption self-energies in 
the n-phonon propagator of the impurity. These can be 
expressed as continued fractions [t|: 
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the latter fraction obviously ending at n-th stage. To 
evaluate the Green's function Gij(u, e) of the full lattice 
problem, we now invert Eq. jSJ with the replacement 

<4 = w - [Hei]i,i — > di=ui - [H e i] izi - T, H<i where 

E H ,i = {GiJ- 1 ~ {Gflr 1 (12) 

is the local self-energy of the Holstein-impurity problem 
at each site i [the off-diagonal terms bi = J + f(ui — Uj+i) 
are kept unchanged]. This treatment of the Holstein 



electron- vibration coupling has formally the same struc- 
ture as Dynamical Mean Field Theory Q, except for 
the absence of a self-consistency loop. It reduces to the 
MA (0) approximation of Refs. [1, |g] at T = and, as 
in MA^ ), it is most appropriate in the non-adiabatic 
regime, Qq > J. 

The above Eqs. (|9l)- (fT2"|) solve the Holstein interaction 
problem for any given configuration of the diagonal and 
off-diagonal disorder. The solution of the full Hamilto- 
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FIG. 1: (a) Spectral function A(k,ui) for a hole in a HOMO 
band, calculated for J = 0.125eV, Xssh = 0.2, too = 0.05 J, 
a% = 0.33, n = 1.5J, 5 = 0.3J, A = 0.2J at T = 0.05J. 
(b) Position of the maxima, tracking the renormalized band 
dispersion ek at low and room temperature (reduced Brillouin 
zone) . 



nian Eq. (JTJ is obtained from Eq. (J7J) upon averaging 
over 50000 realizations of disorder variables on a chain 
of TV = 512 sites. The spectral function is A(k,ui) = 
-iJmG(fc,w), where G(fc,w) = ^ ^ ka{i - j) G itj {oj) 
is the Green's function in momentum space. 

Results. Fig. [TJ, shows a photoemission spectrum cal- 
culated with the above procedure for holes in a one- 
dimensional HOMO band, ideally representing the di- 
rection of maximum conduction in an organic crystal, 
for a representative choice of parameters [slight changes 
in the parameter values within the range indicated af- 
ter Eqs. (l)-(6) do not appreciably modify the scenario]. 
To understand the results it is instructive to analyze the 
different terms in the Hamiltonian separately. We start 
from the interaction with the intra-molecular vibrations, 
H">, that is responsible for the most prominent features 
observed in the spectrum. As is well known, in the case 
of isolated molecules this term gives rise to "shakcofF" 
satellites of the molecular levels, appearing at multiples 
of the vibrational frequency fto [101 ] . The number of vis 
ible satellites is set by the coupling strength, N v n, = a 
and is indicative of the amount of vibrational quanta con- 
stituting the polaronic deformation. 

In a crystalline environment the situation is more com- 
plex. The molecular picture is only recovered when the 
electronic band dispersion is small compared to the vi- 
bration frequency, W <C f2rj. This however does not ap- 
ply to organic crystals, where despite the arguably nar- 
row bandwidths as compared to inorganic semiconduc- 
tors, W is still the largest energy scale in the problem. 
Being W > do, the periodic overtones characteristic of 
molecular spectra are replaced by cuts in the band dis- 
persion. These features are analogous to the "kinks" that 
are commonly observed in the photoemission spectra of 
solids with sizable electron-boson coupling. They add 
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FIG. 2: (a) Temperature dependence of the spectral 
linewidths at different points in the BZ, for the same param- 
eters as in Fig. 1. (b) Widths of the individual subbands and 
total bandwidth. Open circles are obtained including the lat- 
tice thermal expansion calculated in Ref. [24j (the bandwidth 
is normalized to the zero-temperature value of J). 



much affecting the overall bandwidth 0, 0, EH ■ This phe- 
nomenon is clearly visible in Fig. [TJt, where the cosine 
dispersion starting from the top of the band is cut out 
by phonon resonances at Qo and 2f2o (indicated by ar- 
rows). In the case firj ~ (1 — 2) J that is relevant to OSC, 
the interaction with intra-molecular vibrations effectively 
splits the HOMO dispersion into two main subbands sep- 
arated by a sizable direct gap, as shown in Fig. [Tp. The 
gap opens up where the dispersion crosses the first vibra- 
tional cut, which falls accidentally around k ~ 7r/2. This 
result can explain the large separation between the two 
HOMO subbands of pentacene 0, 0] , that is one order of 
magnitude larger than that predicted by ab-initio calcu- 
lations on accounts of the structural non-equivalence in 
the unit cell [HE1- 

Fig. [2k reports the lifetimes for states at different 
points of the Brillouin zone (BZ). The large difference be- 
tween the HOMOl and HOM02 branches, that was also 
observed in Ref. Q , is an additional distinctive feature of 
the intra-molecular interaction ff''': according to the ar- 
guments given above, the HOMOl band lies by construc- 
tion below the threshold for the emission of a vibrational 
quantum, \eu — £fc=o| < ^o, and therefore the electronic 
lifetime there is mostly insensitive to the effects of intra- 
molecular vibrations 0, H| . Vibronic scattering processes 
are instead allowed in the second subband, where they 
cause a much larger line broadening, with linewidths of 
the order of the electronic transfer rate itself. Since only 
the HOMOl states near the band edge can be thermally 
populated or doped in a field-effect device, we argue that 
scattering from high-frequency vibrations should not play 
a predominant role in the transport mechanism of OSC. 

The effects related to if W are essentially temperature 
independent, because the considered vibrations cannot 
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be thermally excited in the relevant range T <C ^o- The 
temperature dependence observed in Fig. [2^, therefore 
arises from other microscopic mechanisms present in the 
model. While the static disorder, , can be excluded 

(it adds a constant broadening oc A to the spectral lines) , 
the scattering from slow inter-molecular lattice vibra- 
tions in if ( M ) is indeed strongly temperature-dependent: 
because the lattice motions are thermally excited for 
T > ojq/2, the linewidths increase linearly with temper- 
ature, which is indeed observed in the HOMOl sector of 
Fig. [5^. Different slopes arise at different points of the 
BZ because the scattering from inter-molecular lattice 
vibrations is minimum at the band edges 

Finally, Fig. 2b illustrates the temperature depen- 
dence of the calculated electronic bandwidths. The total 
HOMO width (black dots) exhibits a moderate increase 
with temperature, that is entirely caused by the coupling 
to low-frequency inter-molecular vibrations: the thermal 
fluctuations of the lattice structure lead to an average 
increase of the transfer integrals given by Jj , = J 2 + a 2 } , 
with aj = AXsshJT for T > uj /2 @, [3, M [2l| . The 
coupling to low-frequency mira-molecular vibrations, not 
considered here, would also lead to an analogous effect 
JUJ . Importantly, this trend is inverted if one takes into 
account the large thermal expansion coefficient charac- 
teristic of organic crystals (open circles in Fig. 2b) (23| : 
upon heating, the lattice expansion leads to an effective 
reduction of the "bare" transfer integral J entering in 
Eq. (2) |24| . The resulting bandwidth reduction quanti- 
tatively agrees with the thermal band narrowing reported 
in pentacene 0, [25[ . This effect should not be confused 
with the the phenomenon of polaronic band narrowing 
[111 , that is the foundation of popular theories of electri- 
cal transport in crystalline OSC, but is inconsistent with 
the parameter range relevant to these materials. 

In standard polaron theories 11 , 2(| 27 1 , that build on 
the molecular limit W <C fio, the dressing of charge carri- 
ers by fast molecular vibrations causes a renormalization 
of the transfer integrals and a consequent reduction of the 
bandwidths, which is exponential in the phonon number 
N v ib- Upon increasing the temperature the "polaronic" 
band further shrinks due to the increased population of 
molecular vibrations, following W — > We~ Nvib ^ T \ with 
N mb {T) = a%[l + 2n B (T)} and n B {T) the Bosc factor. 
For a measurable effect to be observed at room tempera- 
ture, however, one would require fio < lOOO-K" ~ 90meV. 
Such value is considerably smaller than the bandwidth 
and therefore contradicts the molecular limit assumption. 

In the regime of moderate coupling strengths and 
O < W of interest for OSC, it is the very hypothesis 
of an exponential renormalization that fails [H, YL§L ■ 
A polaronic band can still be identified — it is nothing 
but the HOMOl band of Fig.l — but its width is now of 
the order of the vibration frequency fio, and not governed 
by the textbook exponential relation. If we now focus on 
the entire HOMO dispersion, the absence of polaronic 



band narrowing is even more evident: with the present 
parameters one would predict a 30% polaronic reduction 
of the bandwidth from the non-interacting value W = 4 J, 
which is clearly not observed in Fig. 1 (we show in the 
suplementary EPAPS file that this is not a drawback of 
the present method) . The physical reason for the failure 
of the conventional picture is that the molecular degrees 
of freedom are not fast enough to rearrange during the 
electron motion, and therefore cannot provide the instan- 
taneous renormalization implied in the molecular limit. 

Concluding remarks. The calculations presented here 
show that the photoemission spectra of crystalline or- 
ganic semiconductors can differ significantly from the 
non-interacting band picture, reflecting a sizable cou- 
pling of the conduction electrons with other degrees of 
freedom in the system. The most salient features in the 
spectra come from the interaction with high-frequency 
intra-molccular vibrations, that causes the opening of ap- 
parent gaps in the electronic dispersion at multiples of 
the vibrational energy Qq- However, such high-frequency 
modes are unable to efficiently scatter the quasiparticle 
states at the lowest binding energies. For this reason, 
the origin of the temperature dependence of the mobil- 
ity in these materials cannot be explained by polaronic 
effects but should rather be seeked in the interaction 
with inter-molecular vibrational modes of much lower 
frequency, whose temperature-dependent scattering rates 
are directly accessible in ARPES experiments. 

We thank A. Girlando for enlightening discussions. 
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